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~ INTRODUCTION 


This Environmental Science Set of the Hewlett-Packard Computer Curriculum 
Series consists of a set of a Student Lab Book and a corresponding Teacher’s 
Advisor. It was designed to help meet the need for computer-oriented problems in 
environmental science, providing students an opportunity to use the power and 
versatility of the computer to attack some relevant contemporary problems. 





The materials are designed for flexible use as desired by the individual instructor. 
The material in this unit is intended to confront the student with the central 
issues in one of the most pressing problems in the United States today: namely, 
air pollution. The automobile, industry, and overpopulation are discussed as the 
causes of this ecological crisis. The unit emphasizes the ““model”’ approach while 
studying the role of each of these factors in the overall problem. The material can 
be used to supplement and enrich your curriculum in any fashion you choose. For 
example, it can be assigned for individual study or as a class group project. 


The Student Book provides text material, reference data, and programming exer- 
cises for the students. The Teacher’s Advisor contains a brief discussion of the 
important elements of each exercise, as well as sample programs where appropriate. 


Each exercise requires that the students make reasonable assumptions about a given 
pollution problem and then generate a computer model which describes the situ- 
ation. These exercises are geared for general education use for secondary or be- 
ginning college students, but they do vary in degree of difficulty. Some students 
will need a good deal of assistance in order to complete the more difficult prob- 
lems. Though the models are difference equation approximations to systems of 
coupled first order differential equations, the equations are actually very simple to 
solve numerically. Only introductory elements of algebra are required to complete 
the unit and no students should be excluded on the basis of mathematical pre- 
requisites. You will, however, need to see that adequate computer time is available 
to your students. A good planning figure is 30 minutes per students per week. 
Some students will always want to use a great deal more than the allotted time, 
and they should be encouraged if at all possible. 


For best results, you should study all the solutions until you are certain you have 
a complete grasp of the general methods. You will note that there are several 
distinct models presented. You will undoubtedly think of different models. En- 
courage your students to do the same. There are no approved solutions. The more 
initiative the students show in their modelling procedures, the more chance they 
have of arriving at a meaningful solution. 


AIR POLLUTION 
Introduction 
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PROBLEM MODELLING 


The model approach used in this unit is extremely powerful and can be applied 
to a wide range of topics. The student who becomes intrigued with this idea may 
want to push on. An excellent reference is ““The Man-Made World,” Engineering 
Concepts Curriculum Project, Polytechnic Institute of Brooklyn, McGraw Hill 
Book Company. 


Problem Modelling 
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COMPONENTS OF AIR POLLUTION 


The standard nomenclature for air pollution was introduced in this unit, even 
though most of your students have probably not completed a chemistry class. If 
any of your students would like to do some research in general pollution problems, 
you might refer them to Ecology, Pollution, Environment, by Turk, Turk, and 
Wittes, published by W. B. Saunders Co., 1972. 


AIR POLLUTION 
Components of Air Pollution 
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THE AUTOMOBILE AND AIR POLLUTION 


EXERCISE I — Estimating Number of Cars 
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The key to this exercise is that one square mile contains 640 acres. The 
student should estimate the amount of land per dwelling. Remember that 
the area in streets and alleys must be included. Suppose we use the figure 
of 1/2 acre per dwelling. This means that there should be about 1200 
dwellings in the square mile. If we can assume that there are 1.5 cars per 
dwelling, then there should be about 1800 cars in the area. [t doesn't seem 
unreasonable to assume that about 5% or 90 cars are running on the 
average at any time. A good round number to use would be 100 cars 
running atany given time. 


This number is susceptible to measurement if students want to count 
houses and measure traffic flow. A wide range of latitude in assumptions 
still produces something of the order of 100 cars running per square mile 
on the average. Students will more than likely come up with slightly dif- 
ferent numbers, but small differences aren't important. 


EXERCISE 2 — A Simple Model 


Th e algorithni ts 
Pnew = Pold tO RI(IN). 


We begin with P= 0. From Table 2, a car traveling at 40 mph produces 
76.6 ft? of CO per mile. Thus R = (76.6 ft?/mile) (40 miles/hr) or 3064 
ft? /hour. From Exercise 1, N = 100. The students should plot the results 
to see that the algorithm its a linear relationship. 


Program Listing 


100 REM EXERCISE-e2 

110 PRINT 

120 PRINT “TIME",s"'CARBON MONOXIDE" 
130 PRINT "CHOURS)"s"*CCUBIC FEET)" 
140 PRINT 

150 LET P=0 

160 LET R=3064 

170 LET N=100 

186 FOR T=0 TO 24 

198 PRINT TsP 

COO LET P=P+R*N 

216 NEXT T 

999 END 


11 
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Sample Run 

TIME CARBON MONOXIDE 
CHOURS ) (CUBIC FEET) 
QB 0 

1 386400 - 

2 612800. 

3 919200.- 

4 1-22560E +06 
> 1-53208E +06 
6 1-838 40E +06 
7 2-1448GE +06 
8 2-45120E+06 
9 2e7T5760E+06 
10 3eG64QGE+06 
11 3-379 48E +06 
12 3-67680E+06 
13 3-983280E +06 
14 4-28968E+06 
15 4-5960GE +06 
16 4-90240E+06 
17 5 -2688H0E +06 
18 5e¢SIS20E+06 
19 5 -¢82160E +06 
20 6¢12800E +06 
21 6-43448E +06 
22 6-740 80E+06 
23 7-@4720E +806 
24 7¢35360E +86 


EXERCISE 3 — Computing Concentrations 


The volume of air in millions of cubic feet is given by 


If there are P cubic feet of CO then the concentration of carbon monoxide 


V = (5280) (5280) (500)/108. 


in parts per million is given by 


C=P/V 


Program Listing 


100 
112 
120 
139 
142 
156 
160 
170 
186 
190 
206 
210 
228 
999 


REM EXERCISE-3 

PRINT 

PRINT "“TIME'™s"CO CONCENTRATION" 
PRINT "CHOURS)"s" (PPM) 

PRINT 

LET P=0 

LET R=3064 

LET N=100 

LET V=5288*5280*500/1 -O0800E+06 
FOR T=@ TO 24 

PRINT TsP/V 

LET P=P+R*N 

NEXT T 

END 
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Sample Run 


TIME CO CONCENTRATION 
CHOURS) (PPM) 
0 @ 
1 21-9812 
a 43-9624 
3 65-9435 
4 8769247 
2 109-906 
6 131.887 
7 153-868 
8 175-6849 
9 197-831 
10 219.812 
11 241-793 
12 263-774 
13 285-755 
14 307-736 
15 329-718 
16 351-699 
17 373-68 
18 3952661 
19 417-642 
2B 439-624 
21 461-685 
22 483-586 
23 985-567 
24 927 «548 


EXERCISE 4 — Lethal Concentration 


[f the CO concentration is about 500 ppm at the end of one day in our 
model, it should reach 1000 ppm at the end of two days. Obviously, our 
model isn’t accurate, since people would be dying from carbon monoxide 
on a mass basis in residential areas. 


We have assumed that the block of air is 500 feet high. Why not 1000 or 
2000? Also, no air moves in or out in the model, but there must be some 
ventilation. As was pointed out in the introduction to the components of 
air pollution, something unknown ts happening to the carbon monoxide. 
About the only conclusion the student could reach at this point is that 
the simple model can’t be correct! 


EXERCISE 5 — A Garage Problem 


Here the situation is quite different than that in Exercise 4, and much 
more deadly! A production rate of CO of 3064 cubic feet per hour for a 
car operating at 40 mph is equivalent to about I cubic foot per second. 
A reasonable assumption about garage size would be 20 by 15 by 10 feet 
or 3 thousand cubic feet. Thus, if a car engine is operating in our closed 
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garage with an equivalent rate of 40 mph, the CO concentration will pass 
1000 ppm in 3 seconds! 


There most certainly is a hazard in this situation. If the garage ts tightly 
closed with no ventilation, you can be certain that there is a lethal con- 
centration of CO if an automobile engine is running. If you remain in 

the garage long enough for the carbon monoxide to be assimilated into the 
hemoglobin of your blood, that’s it! Depending on the concentration of 
CO, this could happen in a very few minutes. A concentration of 1000 
ppm produces unconsciousness in one hour and death in four hours. How- 
ever, recall that in the garage model, we passed 1000 ppm in 3 seconds! 


This example points out the clear and present danger of any open combus- 
tion taking place in a closed space. Even a charcoal hibachi burning on a 
dining room table can be a threat. Certainly, automobile mufflers that 
leak CO or defective household heating systems that leak CO into the hot 
air system are extremely dangerous. 


EXERCISE 6 — Intersecting Freeways 


We must add the average number of cars expected to be on the freeways 
inside our one square mile residential district. Assume four lanes in each 
direction constitutes a major freeway. Does this seem reasonable? If, on 
the average, the cars are spaced a quarter mile apart (certainly much closer 
during rush hours, and certainly further apart at night), we should expect 
to find 8 cars in each quarter mile section of the freeway, or 32 cars per 
freeway per mile. Since we had two intersecting freeways, we must now 
account for 64 cars operating in the square mile of residential district in 
addition to the 100 in the original model. 


In Exercise 4 we decided that the model was faulty and couldn't be giving 
valid results. However, there certainly is some air pollution problem and 
whatever that is, it is increased by 64% when two freeways intersect in the 
residential district! 


EXERCISE 7 — A Tunnel Problem 
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Let's make the worst case assumptions. Assume there are two lanes of 
traffic moving at 40 mph. At that rate it will take 90 seconds to pass 
through the tunnel. If we assume a car is 15 feet long, and there are two 
car lengths to the car in front, then two cars fone in each lane) are in an 
air space about SO feet long, 25 feet wide (?) and possibly 15 feet high (?). 
This block of air has a volume of about 20 thousand cubic feet. Your 
students will likely make different assumptions but the end results should 
agree reasonably well. 


It makes no difference if we follow the block of air and two cars together 
through the tunnel, or if we examine a fixed block of air and consider the 
cars moving through it. The result is the same in either case. The two cars 
in the block of air are producing CO at the rate of 2 cubic feet per second 
(see Exercise 5). Therefore, in 10 seconds the CO concentration should 
pass the lethal level of 1000 ppm, which means that you would be exposed 
to higher than lethal levels for about 80 seconds if the ventilation system 
on the tunnel went out just as your car entered. 
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Of course, the real hazard takes place if the ventilation system is lost and 
the traffic stalls! [t should be obvious that in such a situation, the motors 
of the stalled cars should be shut off immediately. This exercise should 
stimulate a good discussion in your class. [t certainly points out how you can 
begin to get at seemingly hopeless problems with even a little information. 


EXERCISE 8 — A New Model 


We must make some assumptions to estimate the number of cars in the 
block of airatany given time. If the cars are 15 feet long, separated by two 
car lengths, then about every SO feet on the freeway we will pick up eight 
cars fone in each of the eight lanes). Two thousand feet of freeway would 
hold (2000/50) (8) or about 300 cars. However, we have two freeways so 
there should be about 600 cars in the system. From Table 2, R = 32455 
grams per hour. 


Program Listing 


106 REM EXERCISE 8 

110 PRINT 

120 PRINT “TIME’s"POLLUTANT CONCENTRATION" 
130 PRINT "CHOURS)"s" (MILLIGRAMS PER CUBIC FOOT)" 
140 PRINT 

156 LET P=@ 

160 LET N=60@ 

170 LET R1=3455 

186 LET W=5 

196 LET R2=1-88000E-B2 

200 LET T=2 

210 LET V=2000*2000*500 

220 LET C=1000*P/V 

230 PRINT TsC 

246 LET P=P+R1*N-W*P/50-Re*P 

298 LET T=T+1 

266 GOTO 228 

999 END 
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Sample Run 
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TIME 
CHOURS ) 


omNnoaurnpnwnro®B 


POLLUTANT CONCENTRATION 
(MILLIGRAMS PER CUBIC FOOT) 


8 
1.8365 
1.95898 
2-78 
3e5107 
4.16102 
4-73981 
590625493 
9°71339 
6612141 
6248456 
6688776 
728954 
72635141 
7e¢57925 
7278204 
7296251 
8e12314 
8.26609 
8239332 
8 «50656 
8 «60733 
8 «697803 
877685 
8e84719 
6e91113 
Se 96741 
9eB17TA9I 
9-96207 
9-10174 
9-13785 
9-168 48 
9.19644 
9222133 
9224349 
9226321 
9228075 
9229637 
9.31627 
9232264 
9233365 
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EXERCISE 9 — Equilibrium Concentration 


W 
50 


W 
0 = RN ~ Pog & + R;) 


or Peg = (Ry N)(W/50 + Ro) 


Foy > eg: * Ryan P, POO 


q ed 


| 


and © 1000Peg/V 


Inserting the numerical values of the problem, Cis determined to be 9.36, 
which agrees quite closely with the computer result of Exercise &. 


EXERCISE 10 — Turning the Wind Off 


All we have to do to modify the program from Exercise 8 is to put in a 
counting loop. The new program is shown below. The corresponding 
printout follows. 


Program Listing 


100 
110 
120 
130 
140 
15d 
160 
170 
180 
196 
200 
210 
2208 
23d 
240 
258 
260 
2718 
280 
298 
999 


REM EXERCISE 10 

PRINT 

PRINT “TIME's"POLLUTANT CONCENTRATION" 
PRINT “CHOURS) "s,s" (MILLIGRAMS PER CUBIC FOOT)” 
PRINT 

LET P=0 

LET N=600 

LET R1=3455 

LET W=5 

LET R@=1-9800B0E-92 

LET T= 

LET V=2000*2000*50G0 

LET C=1000*P/V 

PRINT TsC 

LET P=P+R1*N-W*P/50-R2Ox*P 
LET T=T+1 

IF T <= 10 THEN 228 

LET W=0 

IF T >= 21 THEN 999 

GOTO 2208 

END 
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Sample Run 


TIME 
CHOURS ) 


OMNAUHRWN—@B 


Note that at 11 hours the concentrations build up at a steeper rate because 
of the dissipation mechanism (wind) being turned off. If the program is 
continued, a new equilibrium value will be attained. 


POLLUTANT CONCENTRATION 
(MILLIGRAMS PER CUBIC FOOT) 


@ 
120365 
1.95898 
2278 
35107 
4.16102 
4.73981 
9e 25493 
5071339 
6212141 
6-48 456 
6.80776 
7277618 
8e 13492 
9.68407 
10-6237 
11.2554 
12-4749 
13-3867 
14.2893 
15.1829 


EXERCISE 11 — Closing Down Freeways 


We can easily modify the program of Exercise 10 by setting N = 0 instead 
of W = 0 in line 270. Note in the printout that as soon as the source of 
pollutants is gone (N = 0), the dissipation mechanisms start cutting down 


the pollutant concentrations, 





Program Listing 
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188 REM EXERCISE 11 

110 PRINT 

120 PRINT “TIME's "POLLUTANT CONCENTRATION” 
130 PRINT “CHOURS)"*s" (MILLIGRAMS PER CUBIC FOOT)” 
140 PRINT 

156 LET P= 

160 LET N=600 

170 LET R1=3455 

180 LET W=5 

190 LET R2@=1-80800E-82 

200 LET T= 

218 LET V=2000*2000*500 

220 LET C=1000*P/V 

238 PRINT TsC 

240 LET P=P+R1%*N-W*P/50-R2*P 
250 LET T=T+1 

260 IF T <= 10 THEN 220 

276 LET N=0 

280 IF T >= 21 THEN 999 

290 GOTO 226 

999 END 

Sample Run 

TIME POLLUTANT CONCENTRATION 
CHOURS ) (MILLIGRAMS PER CUBIC FOOT) 
Q Q 

1 1.08365 

2 1.95898 

3 278 

4 3e5107 

5 4.16102 

6 4.73981 

7 5 -25493 

8 5.71339 

9 6212141 

10 6248456 

11 6e¢80776 

le 66-8589 

13 5239242 

14 4.79926 

15 4-27134 

16 3-88149 

17 3238333 

18 3281116 

19 2267993 

20 2238514 
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EXERCISE 12 — A Time Dependent Model 


20 


The model itself is easy enough to program. The profile factors must be 
introduced as subscripted variables. This data is read in from DATA state- 
ments in the program. The first number in each DATA statement is X, 
the second Y, and the third Z. 


The students can use any profiles they desire. The set of profiles in the 
program show low traffic during hours of darkness building up to peaks 
in the morning and afternoon. The wind comes up during the late after- 
noon and dies back down in the evening. The pollutant disassociation factor 
is controlled by presence or absence of sunlight. 





106 
110 
120 
130 
149 
15d 
166 
176 
180 
190 
208 
218 
220 
238 
248 
258 
268 
270 
28H 
296 
308 
80 1 
882 
8803 
804 
8805 
806 
807 
808 
809 
812 
S11 
Sle 
813 
814 
815 
816 
817 
818 
819 
826 
821 
822 
823 
824 
999 
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REM EXERCISE 12 

PRINT 

PRINT “TIME's *POLLUTANT CONCENTRATION” 
PRINT "CHOURS)"s"*C(MILLIGRAMS PER CUBIC FOOT)" 
PRINT 

DIM X£L24),YC24)5Z(24) 

FOR I=1 TO 24 

READ XCIJs,YCII)s,Z017] 

NEXT I 

LET P=@ 

LET N=600 

LET R1I=3455 

LET W=15 

LET Re=e05 

LET V=2000*2000*500 

PRINT @5@ 

FOR T=1 TO 24 

LET P=P+R1*XCTI*N-YCT1*W*¥P/5@0-ZCTIJ*RO*P 
LET C=1000*P/V 

PRINT TsC 

NEXT T 

DATA 815051 

DATA «@1 5051 

DATA oB1 sBs 

DATA e@1 5051 

DATA «855051 

DATA elsOaeTl 

DATA o390a0eD 





DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
END 


©75985-3 
©6sG501 
o320a01 
o40B5el 
©4851 
e4selseol 
eAsrelael 
See Sane] 
6929s] 
TS aiaed 
e3ala0e5 
els eSaeT7 
else4s1 
elselsi 
0055051 
O11 5051 
eB1i»Gs1 


2] 
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TIME POLLUTANT CONCENTRATION 
CHOURS ) (MILLIGRAMS PER CUBIC FOOT) 
0 0 

1 ©818365 

2 eeBellTE-Oe 
3 2295662E-B2 
4 3-8 4528E-82 
5 8-e83552E-Be 
6 188913 

7 49514 

8 1.26509 

9 1-88066 

1 2238951 

11 2-19216 

le 3+1928 

13 3249565 

14 3-683083 

15 363841 

16 3225395 

17 3-88633 

18 2234022 

19 1.90396 

26 1.68394 

21 1265287 

ee 1.62205 

23 1655132 

24 1-48 412 


The results clearly show the effects of the various parts of the model. Up 
to about 5 am., nothing happens. Then the morning rush hour and the 
rising sun kicks everything off. Note that we would normally expect to 
get a peak in pollution corresponding to the evening rush hour, but the 
late afternoon wind prevents this. This model is quite powerful and you 
can spend a great deal of time experimenting with your class. [t should be 
particularly useful to make plots varying one parameter at a time to get 
families of concentrations versus time. 


EXERCISE 13 — Political Questions 


to 
tO 


The program from Exercise 14 can be used in this exercise by mrerely 
changing the information in the DATA statements. Have the studeitts 
think a good deal to come up with profile factors that are characteristic 
of the area in which you live. Make reasonable assumptions about numbers 
of cars. It would be informative to put a few additional statements tito 
the problem to permit the program to run for several days of simulation. 
Thus, if all the pollutants are not gone at midnight, the starting value for 
the next day is not zero, and pollutant “stacking” from day to day can be 
observed. Also, students can simulate particular davs which could pose a 
problem, for example, a day without wind. 
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The question of reducing concentrations is designed to force students to 
consider the implications of different courses of action. [t is fine to decree 
that the number of cars will be restricted, or the freeways will be closed 
during pollution alerts, etc., but what are the political and economic ef- 
fects of such decisions? Encourage students to look for realistic decisions 
which improve the pollution problem and which could actually be imple- 
mented from a political and economic point of view. 


23 
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EXERCISE 14 — A Zoning Request 


It certainly would be foolish to approve a zoning request like this. Notice 
that the heavy industry section is surrounded by residential areas and 
parks. We don’t know what type of industries are in the heavy and light 
industry zonings. We don’t know what constitutes a commercial zone. 
There simply isn’t enough information to make any decision at all. 


EXERCISE 15 — Zoning Modification 


The situation is the same as in Exercise 14 except that we are free to 
arrange things to do as well as possible. The point is that this implies some 
criteria as to what constitutes “‘best.’’ The students should decide upon 
reasonable criteria, e.g., keeping the heavy industry as far away from resi- 
dential and park areas as possible, and then use the criteria in locating the 
different functions in the city. 


EXERCISE 16 — More Information 


Students should use their imagination and what has been learned to this 
point in this exercise. What type industries are in the zonings? What 
businesses are in the commercial zone? Is the park and recreation zone 
polluting? What type of pollutants will be expected to come from the 
residential zones? Are there any future plans for annexing to the city 
that could have an effect? 


EXERCISE 17 — A Zoning Model 


First, we must assume typical numbers and types of activities in each of 
the nine segments. Knowing this, we can compute the total production 
of pollutants in each square mile (tons per hour) from the data in Table 3. 
In this solution we will assume: 


LI = Jet Aircraft maintenance + meat packing + 5 electronic plants 
3 tons per day or 0.13 tons per hour. 


R = 100 cars = 760 lbs/hr or 0.38 tons per hour. 

P= = Non-polluting or 0 tons per hour 

HI = Oil refinery + electricity generating plan + railroad shop 
= 64 tons per day or 2.66 tons per hour. 

C = 200 cars or 0.76 tons per hour. 


Your students will surely come up with different assumptions than those 
above. Just check to see that their assumptions are reasonable. 


25 
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Program Listing 


18 

20 

38 

40 

50 

68 

710 

80 

908 

106 
110 
126 
130 
1 48 
158 
160 
170 
180 
196 
200 
210 
220 
230 
246 
259 
260 
270 
280 
298 
300 
318 
320 
3306 
3408 
358 
368 
378 
388 
398 
486 
410 
420 
430 
440 
459 
460 
476 


REM EXERCISE 17 

DIM RC3331,P0353),N0353] 

FOR [=1 TO 3 

FOR J=1 TO 3 

READ RC1IsJ]) 

NEXT J 

NEXT I] 

FOR I=1 TO 3 

FOR J=1 TO 3 

LET PCIs»JI1=8 

NEXT J 

NEXT I 

FOR I=1 TO 3 

FOR J=1 TO 3 

LET NCIsJI=0 

NEXT J 

NEXT I 

FOR T=1 TO 16 

LET NC€i5139=PC015134+RC0151)-.-2*PC1s,1)]+-1*(Pl2,1 J+P 01521) 

LET N€251)=PCe513+RC251)--3*PCeOs1 ]+e-1¥*(PCis,iji+PCles,23+P(351)) 
LET NC351)=P03,114+R(351)]-.-2*P(351)+-1*(P£2513+P035,2)) 

LET NCis2)=P01,2)4+R015,2)--3*PC1,20)+-1*(P£151)1+PCes,2e+Plls3])) 
LET NC2,2)3=P0252)+R602521)--.4*P02,2)+-1* (P0251 +P Cis2)+Pl2,31+P03,2e)) 
LET NC352)=P(03,2)+R(352)-.3*P(3,2)4+.1*(P 0351 1+PC2,20)+P0353)) 
LET NC153)=P0153)4+R01531-.2*P0153)+-1*(P01,23+PC0253)) 

LET NC253)=P02531)4+R0253)--3*P02,3)+.1* CP L2,2I14+PC1s»3I+P0353)) 
LET NO3533=P0(353)4+RO03,33)--24P(353)+-1*(P03,2)+P C2531) 

FOR I=1 TO 3 

FOR J=! TO 3 

LET PCI,JIJ=NCI,J] 

NEXT J 

NEXT I 

PRINT 

PRINT "FT = “sT 

PRINT 

FOR I=1 TO 3 

FOR J=1 TO 3 

PRINT PCI,J]); 

NEXT J 

PRINT 

PRINT 

NEXT I 

NEXT T 

DATA «133+-385-76 

DATA 032-665-38 

DATA 2385-38513 

END 





Sample Run 


T= 1 


e13 


Q 


38 


e2te 


0317 


©7122 


T = 8 


©4794 


©7583 


1.0856 


T= 4 


75828 


1.24833 


1.48374 


T= 5 


1.09976 


1-75799 


1691413 


T= 6 


1.4919 


2227646 


2e370713 


1.6893 


5°6182e 


125943 


2238301 


659962 


2222383 


3286295 


1¢44473 


2283628 


3-72554 


6-280131 


3243393 


1.444 


1-001 


o3l 


2e1154 


11-6931 


S744 


20719056 


2039517 


91826 


3-47027 


3 «88746 


1.32643 


4.15126 


376537 


1678352 


AIR POLLUTION 
Industry and Air Pollution 
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T= 7 

1.92372 4237233 4.8301 
2019991 Be 9IBBIC 4242937 
2<«847162 4.01931 2e2Tl6T4 
T= 8 

223862 9+8061 9058425 
3e32717 9-56264 3298133 
3-34902 4.59604 279626 
T= 9 

2087229 5262958 6217214 
368579 10-1986 5e72325 
3e84434 5216712 3¢33475 
T = 10 

3237658 6224561 6-833 
4.39206 10.817 6 «35683 
4-35797 5°73476 388683 





The matrix R has pollutant production rates in each segment correspond- 
ing to the activities as outlined above. Matrix P is the old level of pollutants, 
and matrix N is the new level. Lines 160 through 240 apply the algorithm 
to each of the nine segments. The balance of the program follows with- 
out difficulty. 


The sample run shows that the pollutants from the heavy industry seg- 
ment diffuse out into the other segments as the hours increase. Since there 
are no losses in the system as a whole, the pollutants will keep on increas- 
ing in each segment and will gradually approach the same concentration 
at the same time in each of the segments. At this point, the city is “full.” 
Of course, losses should be built in to make the model more reasonable. 
You could easily assume a wind loss adjustment to the model as was done 
previously. Encourage your better students to experiment with the model. 





AIR POLLUTION 
Industry and Air Pollution 
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EXERCISE 18 — Restructure of the City 


This is an open-ended exercise. Students must come up with some criteria 
to use in judging the worth of a given distribution of activities. An example 
of such a criteria might be to keep the pollutant levels as low as possible 
for as long as possible in the residential areas. With such an assumption, 
students can shift the activities around by modifying the DATA State- 
ments in the program of Exercise 18, running the program and analyzing 
the results. 


EXERCISE 19 — Model Your Own Town 


This is another open-ended exercise. Here your students must apply most 
of the knowledge they have acquired in preceding exercises. You should 
make sure that assumptions are reasonable and that the resultant model 
corresponds fairly well to the real town that is being modelled. [t may not 
be possible to make a rectangular model of your town in one-mile squares, 
so feel free to modify the basic model. 
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PEOPLE AND AIR POLLUTION 


EXERCISE 20 — People and Pollution 


Since we have no dissipation mechanisms, we may as well consider the 
pollutant concentration per square mile. We will consider the air over each 
square mile up to an altitude of 500 feet. The program follows without 


difficulty. 


Program Listing 


Sample 


180 

112 

128 

13@ PRINT 
140 PRINT 
150 PRINT 
162 

170 PRINT 
180 

199 LET 
200 LET 
21@ 

220 

230 

240 NEXT N 
999 END 
Run 


POPULATION = 


PART 


5° 86268E-94 


POPULATION = 


PART 


1-1 7252E-33 


POPULATION = 


PART 


1-75878E-33 


POPULATION = 


PART 


2-34504E-33 


POPULATION = 


PART 


2-9313BE-83 


"POPULATION = 


REM EXERCISE 26 
LET V=5280*5280*500 
FOR N=S880 TO 25000 STEP 5800 


5888 DAILY 
ORG 


6-45661E-92 


18000 DAILY 
ORG 
129132 
15808 DAILY 
ORG 
193698 
2O0BGB DAILY 
ORG 
2258264 
25808 DAILY 
ORG 
322831 


PRINT "PART'’,"ORG''s’*SO''s*NO''s "CO" 


LET P1=(N/25)*9 -O008B0E-02*454800-/V 
01=(N/25)#*4-5/(V/1 -OOBGBE+O6) 
S1=(N/25)*-.7/CV/1 -BBDOOGE +06) 
LET N1=(€N/25)*2-6/0V/1 -OBO00B0E+86) 
LET C1=(N/25)*32.6/(V/1 -@BBB0E +86) 
PRINT P1,01551,N1,C1 


RATES PER SQUARE 


18) 


1-930436E-B2 


RATES PER SQUARE 


SO 


2-3887T2E-B2 


RATES PER SQUARE 


SO 


3-8@1389E-82 


RATES PER SQUARE 


50 


RATES PER SQUARE 


SO 


5-92181E-82 


"S3N3"DAILY RATES PER SQUARE MILE" 


MILE 


NO 


3-TIG4IE-B2 


MILE 


NO 


7+46097TE-82 


MILE 


NO 


0111915 


MILE 


NO 


°149219 


MILE 


NO 


186524 





CO 


«467746 


CO 


935491 


CO 


1.40324 


CO 


1.87898 


CO 


2.33873 


People and Air Pollution 
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EXERCISE 21 — Dissipation Mechanisms 


Open-ended problem — turn your students loose! 


EXERCISE 22 — A Segmented Model 


Open-ended. Use some of the approaches developed previously. Encourage 
imagination and hard work. 


EXERCISE 23 — A Real City 


The point here is that, as one proceeds inland from the ocean, the air 
pollution problem worsens and approaches some kind of bad limit as you 
approach Riverside. If the model is reasonable, it must account for this! 


a2 
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THE FUTURE 


EXERCISE 24 — Pollution In Your Home Town 


This is a good exercise in civics and air pollution. Students should be very 
curious and fairly knowledgeable about air pollution at this point. 


EXERCISE 25 — Pollution Control Programs 


This exercise follows from Exercise 24. 


The Future 
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